Do we know the mass of a black hole? Mass of some cosmological black hole models 
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Using a cosmological black hole model proposed recently, we have calculated the quasi-local mass 
of a collapsing structure within a cosmological setting due to different definitions put forward in 
the last decades to see how similar or different they are. It has been shown that the mass within 
the horizon follows the familiar Brown- York behavior. It increases, however, outside the horizon 
again after a short decrease, in contrast to the Schwarzschild case. Further away, near the void, 
outside the collapsed region, and where the density reaches the background minimum, all the mass 
definitions roughly coincide. They differ, however, substantially far from it. Generically, we are 
faced with three different Brown- York mass maxima: near the horizon, around the void between the 
overdensity region and the background, and another at cosmological distances corresponding to the 
cosmological horizon. While the latter two maxima are always present, the horizon mass maxima 
is absent before the onset of the central singularity. 

PACS numbers: 95.30.Sf,98.80.-k, 98.62. Js, 98.65.-r 
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I. INTRODUCTION 



What does general relativity tell us about the mass 
of a cosmological structure in a dynamical setting? We 
know that massive sources produce gravitational field 
which has energy. In relativity theory, the equivalence 
of mass and energy means that it is only the combined 
energy which may be measured at a distance. We should, 
therefore, expect that because of non-linearity of the 
gravitational field the mass of the material source, its 
kinetic energy, and the gravitational energy it produces 
combine in a nonlinear and non-local way to produce the 
effective energy. In simplest case of spherical symmetry 
in vacuum, this effective energy is just the Schwarzschild 
energy. Although this argument about the effective 
energy seems very plausible, there are still disputes in 
the literature simply about the definition of mass in the 
spherically symmetric vacuum cases i(for an interesting 
discussion about the concept of mass in relativity see 
also 0). What if the material mass is embedded in a 
cosmological setting? 

In a cosmological setting, due to the asymptotically 
non-flatness of the cosmological background, one can not 
use global definitions such as ADM mass and angular 
momentum [lj. This has led authors in the last decades 
to the notion of the quasi local mass (QLM) or quasi 
local energy (QLE), applicable to non-local structures 
in any general dynamical situations 2] . In this paper 
we will interchangeably use QLM or QLE as the same 
concept, although they may differ in other contexts. It 
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is a fundamental fact of general relativity reflected in 
the equivalence principle that there is no such concept 
as the local mass of a local object: a local object at 
the origin of a freely falling frame will not experience 
any gravitational acceleration. In the absence of a local 
gravitational effect the notion of a QLM, corresponding 
to a two dimensional compact surface, although not 
uniquely defined, is the only concept one may try to 
apply to astrophysical cases. 

Let us define a cosmological structure as an overdensity 
region within a cosmological background which we 
assume to be asymptotically FRW. There are not much 
viable exact models representing such a structure. We 
will rely on a recent analytical model proposed by 
the authors based on a inhomogeneous cosmological 
LTB model to construct such an asymptotically FRW 
universe including an overdensity region evolving to a 
black hole different from that of the Schwarzschild Q. 
Our aim is to understand the notion of mass due to 
different definitions of QLM for such a cosmological 
structure within a FRW universe. 



In general, there is no unique definition of the quasi 
local energy, except at infinity in an asymptotically flat 
space-time, where one has the Arnowitt-Deser-Misner [l[ 
energy Eadm at the spatial infinity and the Bondi-Sachs 
§ energy Ebs at the null infinity. Hawking [l2| defined 
a quasi local mass which has various desirable properties: 
it is zero for a metric sphere in flat spacetime, gives the 
correct mass for the Schwarzschild solution on a metric 
sphere, and tends to the Bondi mass asymptotically in a 
static asymptotically flat spacetime. However, it is non- 
zero for generic 2-surfaces in flat spacetime. Hayward [l4| 
has proposed an expression for quasi local energy which 
maybe considered as a modification of Hawking's energy 
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by shear and twist terms. Another attempt was due to 
Misner and Sharp [|| with a well-understood Newtonian 
limit OH. 

The most promising QLM definition, however, seems to 
be the one proposed by Brown- York [l5|. Motivated by 
the Hamiltonian formulation of general relativity, they 
found an interesting local quantity from which the defi- 
nition of quasi local mass was extracted. Their definition 
depends, however, on the choice of the gauge along the 
3-dimensional spacelike slice. It has the right asymptotic 
behavior but is not positive in general. Motivated by 
some geometric consideration, Liu-Yau 11811 (see also Ki- 
jowski Booth-Mann (l6| . and Epp [l7| ) introduced a 
mass which is gauge independent, and always positive. 
However, it was pointed out by O'Murchadha et al [13] 
that the Liu-Yau mass can be strictly positive even when 
the 2-surface is in a flat spacetime. 

Our aim is to study some of these quasi local masses for 
spherically symmetric structures in a cosmological set- 
ting within general relativity to see how similar or dif- 
ferent they are. We review different definitions of quasi 
local masses in section II, followed by the introduction 
of LTB metrics and the corresponding mass definitions 
in section III. In section IV we explicitly calculate nu- 
merically the masses for two cosmological toy black holes 
and a structure with an NFW mass profile [23[ . We then 
conclude in section V discussing the results. 



II. QUASI LOCAL MASS DEFINITIONS 



assuming the energy momentum tensor for the perfect 
fluid in the form 

T* = -p(t,r), T;=p r (i,r), T e e = T$ = Pe (t,r), (2) 

with the week energy condition 

P>0, p + Pr >0, p + Pe >0, (3) 

we then obtain the Einstein equations in the form 



2M' 



P = 



R 2 R 



7 , Pr = 



2M 

rzr' 



2{p e - p r ) R' p' r 



p + Pr R P + Pr' 

-2R' + R'^ + R^ = 0, 

where 

G = e- 2 ^{R') 2 , H = e- 2u {R) 2 , 
and M is defined by 

2M 



G-H = 1 



R 



(4) 



(5) 



(G) 



(7) 



(8) 



The function M can also be written as 



In general relativity the mathematical entity used to 
define the mass-energy is the symmetric energy momen- 
tum tensor T^ v , representing the source-current for grav- 
ity, although the proper interpretation of it is only the 
source density for gravity. This fact is one of the roots of 
difficulties to define mass in general relativity. Reason- 
able total energy-momentum can be associated with the 
whole space-time provided it is asymptotically flat. This 
has led general relativist to 'quasi-localization' of total 
quantities, and construction of 'quasi-local' mass-energy. 
Techniques used in the quasi localization depend on the 
actual form of the total quantities, yielding inequivalent 
definitions for the quasi- local masses Q. Here we out- 
line some of the mostly used definitions before going on 
to apply them to specific models for mass condensation 
within FRW cosmological models and try to interpret the 
results. 



A. Misner-Sharp mass 

Take a collapsing ideal fluid within a compact spher- 
ically symmetric spacetime region described by the fol- 
lowing metric in the comoving coordinates (t,r,9,if): 



M = - I pR dR, 



M =T I P 
8t Jo 



where 



and 



d 3 V = Ane^R'dr, 



d_ 

d^ 



dt 



(9) 



(10) 



(11) 



(12) 



ds l 



(1) 



The last form of the function M indicates that when 
considered as energy, it includes contribution from the 
kinetic energy and the gravitational potential energy. M 
is called the Misner-Sharp energy. 

Hayward [To| showed that in the Newtonian limit of a 
perfect fluid, M yields the Newtonian mass to the lead- 
ing order and the Newtonian kinetic and potential en- 
ergy to the next order. In vacuum, M reduces to the 
Schwarzschild energy. At null and spatial infinity, M 
reduces to the Bondi-Sachs and Arnowitt-Deser-Misner 
energies respectively [loj . 
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B. Hawking mass 

Hawking |l2| defined a quasi-local mass for the space- 
like topological 2-sphere S: 



E„(S) = 



Area(S) 1 
16ttG 2 1 2tt 



(*(-*2-trA + $ii+A)d5), (13) 



^reo(S) 
16ttG 2 



where the spin coefficients /? and p' measure the expan- 
sion of outgoing and ingoing light cones. For the defi- 
nitions of o") A, $n and A see [Hj]. The Hawking 
mass has various desirable properties: it is zero for a 
sphere in flat spacetime, gives the correct mass for the 
Schwarzschild solution on a metric sphere, and tends to 
the Bondi mass asymptotically in a static asymptotically 
flat spacetime. It is invariant under the boost gauge 
transformation. It is, however, non-zero for generic 2- 
surfaces in a flat spacetime (see for more detail). 



C. Hay ward mass 



Hayward, using a 2+2 formulation of general relativity, 
jives the following definition for a quasi- local energy [14j : 



term that is used to normalize the energy with respect 
to a reference spacetime, not necessarily flat. This quasi 
local mass has the right asymptotic behavior but is not 
positive in general. The BY mass depends, however, on 
the choice of the gauge along the 3-dimensional spacelike 
slice £. To avoid this gauge dependence , Yau [18[ (see 
also [HI] and [l7j ) introduced a mass which is gauge inde- 
pendent and always positive. Based on Yau's definition, 
Liu- Yau defined the mass 



E = — 



1 

8^ 



(fx^(s/k 2 -I 2 - k ), 



(16) 



where I and k are traces of extrinsic curvatures l a b = 
cr a a fflc u <i an d k a b = a£afV c n<j respectively, for a ab = 
9ab + u a ui, — n a nb being the metric on the 2-sphere B. 
It was pointed out [20l | that this Liu- Yau mass is strictly 
positive, even when the surface is in a fiat spacetime. 
It can be shown that for static spherically symmetric 
spacetimes the Brown- York quasilocal energy at the sin- 
gularity is zero 22] , in contrast to to the Newtonian grav- 
ity, in which the energy of the gravitational field diverges 
at the center for a point particle. Apparently, the nonlin- 
earity of general relativity has removed this infinity. This 
QLE attains its maximum inside the horizon having an 
infinite derivative just on the horizon, before matching 
to its value outside the horizon: the black hole looks like 
an extended object. 



E H (S) 



Area{S) 1 
16ttG 2 [ 2tt 



-O a bO " 



-u a u a )dS\, 
(14) 



where a ab and a ab are shears and oj a is the normal fun- 
damental form. The energy is zero for any surface in 
flat spacetime, and reduces to the Hawking mass in the 
absence of shear and twist. For asymptotically flat space- 
times, the energy tends to the Bondi mass at null infin- 
ity and to the ADM mass at spatial infinity. It depends, 
however, implicitly on the gauge choice [2j. 



D. Brown- York mass 

Looking into the Hamiltonian formulation of general 
relativity, Brown and York [15| found interesting local 
quantities from which the definition of quasilocal mass 
was extracted. Consider a 3-dimensional spacelike slice 
£ bounded by a two-surface B in a spacetime region that 
can be decomposed as a product of a spatial three-surface 
and a real line-interval representing the time. The time 
evolution of the two-surface boundary B is the timelike 
three-surface boundary 3 B. When £ is taken to intersect 
3 B orthogonally, the Brown- York (BY) quasilocal energy 
is defined as: 



E = — 
8tt 



— / d 2 xy/a{k 
i7r Jb 



ko), 



(15) 



where a is the determinant of the 2-metric on B, k is the 
trace of the extrinsic curvature of B, and fco is a reference 



III. LTB METRIC 

The LTB metric may be written in synchronous coor- 
dinates as 

ds 2 = dt 2 Trrdr 2 - R(t,r) 2 dn 2 . (17) 

1 + f[r) 

It represents a pressure-less perfect fluid satisfying 

Here dot and prime denote partial derivatives with re- 
spect to the parameters t and r, respectively. The an- 
gular distance R, depending on the value of /, is given 
by 

M 

fl=-y(l-cosfa(r,t))), 
V-sm(r,) = ^—(t-t n (r)), (19) 



for / < 0, and 



for / = 0, and 



AI 



R=(-M)*(t-t n y- 



M 



(20) 



r =— (cosher, t))-l), 
smh{r,)-n= f -^{t-t n {r)), (21) 
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for / > 0. 

The metric is covariant under the rescaling r — > f(r). 
Therefore, one can fix one of the three free parameters of 
the metric, i.e. t n (r), f(r), and M(r). 
This metric has two generic singularities: the shell fo- 
cusing singularity at R(t, r) — 0, and the shell cross- 
ing one at R'{t,r) = 0. However, if ^f RI and are 
finite at R = then there is no shell focusing singu- 
larity. Similarly, if ^ is finite at R' = Q then there 
is no shell crossing singularity. To get rid of the com- 
plexity of the shell focusing singularity, corresponding to 
a non-simultaneous big bang singularity, we may assume 
tn{r) = 0, which will be the case for our toy models. This 
will enable us to concentrate on the behaviour of the col- 
lapse of an overdensity region in an expanding universe 
without interfering with the complexity of the inherent 
bang singularity of the metric @. 



ft is easy to show that 



>(£)\R=2M 



0. Therefore, there 



may exist an apparent horizon, defined by R = 2M, be- 
ing obviously a marginally trapped tube. It will turn out 
that this apparent horizon is not always spacelike and 
can have a complicated behaviour for different r [6J. 



A. Misner-Sharp mass 

It is easily seen from © that M(r) in the LTB metric 
is identical to the Misner-Sharp mass. The rate of change 
of this mass for any R = const in the collapsing region is 
positive as can be seen by the following argumentation. 
Noting that R < in the collapsing region, and assuming 
no shell crossing, R' > 0, we obtain from R'dr + Rdt = 



that &Ih= 



di 



\R— const 



> 0. Therefore, given 



dM(r) 
dr 



> 0, we see 



, i , dM (r) I _ dM (r) dr I n 

tnat — — — \R=const — — -JZ — -JJ\R=const > U. 



dt 



C. Brown- York mass 

The 2-boundary B maybe specified by r — constant 
and t — constant. We then obtain for the trace k o f the 
extrinsic curvature fc Q f, for LTB's metric k = — 2 ^ + ^ ■ 



The Brown- York energy is then given by 

Mby = —R\/l + / — Subtraction term. 



(25) 



The subtraction term is chosen to be the correspond- 
ing FRW term for the t = constant slice, i.e. 
—R\/l + J\frw- Now the rate of change of the Brown- 
York mass is given by 



dM 



BY 



dt 



-Ry/T+f + (R^T+f)\ 



FRW- 



(26) 



The first term is responsible for the flow of dust falling 
into the center and the second term is due to the cosmo- 
logical expansion. In the collapsing phase of the central 
region the first term is dominant and the Brown- York 
mass increases within the sphere of radius r. 
The Brown- York mass may also be calculated for the 
2-boundary B with the constant physical radius R — 
constant at t = constant, leading to 



2M 

MBy\R.=const — 1 Subtraction term. 

(27) 

The subtraction term is again the corresponding FRW 
term as the background. 

Similarly, the Liu-Yau mass for the LTB black hole model 
is given by 



2M 

Mly = 1 Subtraction term, (28) 



B. Hawking and Hay ward masses 

The LTB null tetrad needed to calculate the Hawking 
mass is given by 



VTTf 



^7r^-<>.<>i. n" = (i, -V±±1,0,0), (22) 



and 

- ^(0,0,1,^-), m^ = -^(0,0,1,^). 
i?V2 v szn6» ; ' R^/2 K sine' 

(23) 

We then obtain for the Hawking mass 

M Haw = M(r), (24) 

i.e. it is equivalent to the Misner-Sharp mass, as ex- 
pected. This is due to the vanishing of twist and shear in 
metrics being spherically symmetric for round 2-sphere. 
We also conclude that the Hayward mass is identical to 
Hawking mass for the LTB metrics. 



which is valid for k > I, or for the untrapped region 
R > 2m and Subtraction term=— R . The rate of change 
of this mass is given by 



dt 



dM LY ■ I 2M RM d< Subtraction term) 
-R\ I 



R 



dt 



(29) 

which is again an increasing function within the collaps- 
ing region. 



IV. MASS OF EVOLVING BLACK HOLES 
WITHIN FRW UNIVERSE 

We are now interested in the mass of a cosmological 
overdensity region evolving into a black hole in a FRW 
background. We first choose two toy models and look for 
the mass of structures they represent. This should give us 
an overall view of the different mass definitions and their 
differences. Then we go to a more realistic model start- 
ing with a given density profile and look for the model 
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parameters and the corresponding masses. Our cosmo- 
logical black hole is going to be modeled by a LTB so- 
lution representing a collapsing overdensity region at the 
center and a flat FRW far from the overdensity region Q . 
The overdensity region may take part in the expansion of 
the universe at early times but gradually the expansion 
is reversed and the collapsing phase starts. For a more 
realistic model we assume the familiar NFW profile [23[ 
for the overdensity region within a LTB model and look 
for its consequences as regards different mass definitions. 



A. Example I:Toy model I with lim^oo f(r) — > CP; 
structure within an asymptotically closed-flat LTB 
metric 

The model is defined by the requirement f(r) < and 
f(r) -> when r -> oo, and M(0) = 0. In both LTB 
toy models we assume t b = 0. Let us use the ansatz 
f(r) = — re~ r leading to 



M{r) 



r 3/2 



(1 



r 3/2 



where a is a constant having the dimension [a] = [L]~ 2 
@. The constant a is fixed by ato = 3tt/2, corresponding 
to the collapsing mass condensation around r = start- 
ing in the expanding phase of the bound LTB model. 
Equations ([TO)) and ([20)1 then lead to 



(30) 



R= V^(l + r 3 / 2 ) 



(f - cosrj(r,t)), 



n — sin(rj) 



(l + r 3/2) 



at. 



Fig-© shows schematically the behavior of the curvature 
function f(r) and the corresponding Brown- York mass 
for a sphere of constant co-moving radius r = constant. 



R = 



r 3/2 



a(e 



-(1 — cos77(r, t)), 



j] — sinrj - 
for / < and 



r 2 +20000' 



r 2 +20000 



1.5 



(l + r 3/2) 



-at, 



(31) 



R= , V \ i-r(cosh»7(r,t)-l), 



r 2 +20000 



I ± p-rU.5 

• > V r 2 + 20000 I , 

* ~ SmhT > = (l + r 3/2) Q< ' 



(32) 



for / > 0. 

Fig-© shows schematically the two Brown- York 
masses inside spheres of constant comoving radius r for 
the closed-flat and open-flat cosmological black hole toy 
models. Note the negative values of the BY mass in the 
case of open-flat model at distances far from the central 
overdensity region. In the case of closed-flat model, the 
Brown- York mass behaves similar to the corresponding 
Schwarzschild mass 22]. In both cases the Brown- York 
mass is zero at the central singularity in contrast to the 
Misner-Sharp and Hawking and Hayward mass, and re- 
mains finite within and on the horizon. 



f(r) 



ll-FRW: k=0 



l-FRW: k=0 



B. Example II: Toy model II with lim r ^oo f(r) — >• + ; 
structure within an asymptotically open-flat LTB 
metric 

What would happen if we choose the curvature func- 
tion f(r) such that it is negative for small r but tends 
to zero for large r while it is positive? We still have 
a model which tends to a flat FRW at large distances 
from the center, having a density less than the critical 
one corresponding to an open FRW model. The model is 
defined by the ansatz f(r) = —r(e~ r — r „ 1 +c ) with n = 2 
and c = 20000, leading to 



M(r) 



1 



.3/2 



(1 



r 3/2 



Asymptotically flat model 1 
Asymptotically flat model 2 



FIG. 1: The upper diagram shows curvature functions f(r) 
for the closed-flat (I) and open-flat models (II). The lower di- 
agram shows schematically two Brown- York masses for con- 
stant comoving radius r. 



where a is a constant having the dimension [a] = [L]~ 2 
and fixed by the requirement at = 37r/2. Equations (|19[) 
and (|20p then lead to 



The more interesting case of the mass inside a sphere 
of constant physical radius R is shown in Fig. ^ just for 
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the closed-flat case. The corresponding mass behaviour 
for the open-flat case is similar to the closed-flat one, 
at least in the range of radius plotted in Fig.Q. Note 
the Q— value (f2 = p/p cr ) as a function of R plotted 
in the same figure. A comparison to the BY mass of 
the Schwarzschild metric [22| shows following similarities 
and differences. The Brown- York mass of our cosmolog- 
ical black hole toy model has two maxima: one near the 
apparent horizon and the other far from the horizon at 
cosmological distances, corresponding to the FRW cos- 
mological horizon, and at the same time around a void 
occurring in the models studied. This void is, however, 
too shallow to be seen in the figure. We will see in the 
next section that in more realistic cases this second max- 
imum splits in two different ones: one around the void 
where the density tends to the background value 51 = 1, 
and the other around the cosmological horizon which is 
out of the range of our interest. The first maximum oc- 
curs just before the apparent horizon of the central black 
hole followed by an infinite slope similar to the BY mass 
of the Schwarzschild metric [22| . Let us call it the Hori- 
zon mass maximum, which can be seen to be related to 
the non-zero Misner- Sharp mass of the black hole. After 
reaching a minimum, the BY mass then increases again 
with increasing physical radius R up to the second max- 
imum, in contrast to the BY mass for the Schwarzschild 
case which decreases to the ADM mass. This second 
maximum, however, occurs after an infinite slope just af- 
ter the cosmological particle horizon, in contrast to the 
behavior of the horizon mass maximum which has the 
infinite slope before the mass maximum and the appar- 
ent horizon. Because of the fact that this maximum oc- 
curs around the region where the density have reached 
the background FRW density, after passing a void and 
separating the central black hole from the almost homo- 
geneous background, we call it the structure mass max- 
imum. In contrast to the one corresponding to cosmo- 
logical distances which we may call cosmological mass 
maximum, the separation of these two maxima will be 
obvious in the example III. Note that The Liu-Yau mass 
behaves similar to the Brown- York mass in the strong 
gravity region, but it approaches the Misner-Sharp mass 
far from the center. 

Fig-© shows the BY mass at three different times, just 
before the onset of singularity and after the singularity 
has appeared, corresponding to different density profiles. 
It is obvious from the figure that the mass maxima in- 
crease and shift towards larger R values as the time in- 
creases. The Horizon maximum is missing before the 
onset of the singularity. This maximum may be used as 
an indicator of singularity appearance in numerical rela- 
tivity. The structure mass maximum is, however, present 
for any density profile irrespective of the occurrence of a 
central singularity (see Fig flU)) . 



— -MiSh BrYo-r BrYo-R — LiVa Q 




FIG. 2: Different masses of the cosmological black hole toy 
model I for constant physical radius R. Horizon radius is 
about 0.1 in terms of 7? units. 




• BrYo-l — BrYo-ll BrYo-lll G-l 



FIG. 3: Brown- York mass for three density profiles corre- 
sponding to three different times after the onset of the sin- 
gularity and a time before the singularity has been formed. 
As the mass in-fall increases the black hole mass maximum 
increases too and shifts towards larger R values. For dashed 
line, horizon radius is about 0.1 in terms of R units. 



C. Example III: Structures within asymptotically 
flat LTB models having a NFW density profile 

For a more realistic modeling of a cosmological struc- 
ture, we need an algorithm to construct functions M, f, 
and tb in the LTB solution from physical quantities of 
the system such as the density profile. Krasihsky and 
Hellaby [24| propose an algorithm by which knowing the 
initial and final density profiles of an object one can find 
/ and tb as functions of M. For the sake of simplicity, it 
has been assumes here M = r. 

We choose a Gaussian profile for the density at the last 
scattering surface as the initial time. The final profile, 
say at z ~ 0.2, is then chosen to be the universal density 
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profile for the dark matter suggested by Navarro, Frenk 
and Wright [2J|. To simulate a void compensating the 
overdensity mass region relative to the cosmological back- 
ground [25j | r .'2( j , we convolute density profile of structure 
with a Gaussian underdensity. At far distances from the 
center of structure, density tends to the critical density, 
corresponding to a flat matter dominated cosmological 
background: 



p;(r) 

PNFw(r) 
Pf( r ) 



I T \2 

Pcrit(U)((ScMBe kR " ; 
Sc 

Pcrit ~ 



Pcrit (*/)(( 



:) 2 

Sr. 



b 2 )e 



1) 



1) 



Using this algorithm, we have calculated the correspond- 
ing LTB functions, needed to define the metric and to 
calculate different masses. / as a function of M is de- 
picted in Fig. (0| . Note that, in contrast to toy models 
discussed above, now the bang time tb is non-vanishing 
as is shown Fig.©. Theses functions have been calcu- 
lated for r-values larger than a minimum corresponding 
to the central singularity R — 0. To calculate the mass 
we have assumed r s s» 500fcpc and S c « 4, 000 [27]. 




1E+11 1E+13 1E+15 1E+17 1E+19 1E+21 1E+23 1E+25 1E+27 
M(M ) 

FIG. 4: f(M) for a cluster of galaxies with NFW density 
profile in a fiat background. 



The result for different mass definitions is depicted in 
Fig.© as a function of the physical distance from the 
center. The mass up to the void turns out to be of the 
order of 10 14 solar masses, which turns out interestingly 
to be of the same order for all mass definitions. 

There are some interesting features in the Fig.©. 
Note first that the horizon mass maximum, occurring also 
for the NFW profile, is not depicted in the Fig.® as it 
is very near to the origin of the figure. Numerical calcu- 
lation shows that the apparent horizon is at a physical 
distance of the order of 10 n fcm = O.Olfcpc correspond- 
ing to a Misner-Sharp mass of the order of 3 x 1O 11 M . 
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1E+11 1E+13 1E+15 1E+17 1E+19 1E+21 1E+23 1E+25 1E+27 

m(m ) 

FIG. 5: tt(M) for a cluster of galaxies with NFW density 
profile in a fiat background. 



The structure maximum is shifted more to the left at 
distances of the order of few Mpc corresponding to the 
place of the void just before the density profile reaches 
the background values. Different mass definitions seems 
to coincide at this structure maximum which is about 
10 14 for the model we have constructed and it is about 
300 times the horizon maximum mass and about 30 times 
the mass up to the distance of about IMpc, correspond- 
ing to CI « 100. 

The cosmological mass maximum of Brown- York, which 
appears at cosmological distances, is now separated from 
the structure mass maximum. It is interesting to note 
that BY mass for constant co-moving and physical ra- 
dius remains almost the same for distances after the void, 
while the Misner-Sharp and Liu-Yau masses increase, be- 
ing almost equal. Fig.© shows the corresponding dia- 
gram for galactic masses. We recognize similar features 
as those for cluster masses, except the less exposed void 
at distances less than one Mpc. 

The rate of matter flux through the apparent horizon 
(as a quasi-local black hole boundary in an asymptot- 
ically flat universe) is a useful quantity in astrophysi- 
cal studies of black holes as well as theoretical study of 
black hole laws [28|. The matter flux for Misner-Sharp 



dM(r) 
dt 



AH 



mass along the apparent horizon is given by 

dJ rfr 1&AH' ^ s we see * n Fi§-®' the ra t e °f the matter 
flux increases with time at the initial phase of the black 
hole formation, up to a maximum value of the order of 
magnitude 10 6 solar masses. It then decreases while the 
collapse is continuing and the black hole boundary (ap- 
parent horizon) freezes out in the expanding background 
(see Fig ©. In the case of BY mass, the rate of matter 
flux is similar to that of Misner-Sharp one (see Fig.©), 
being almost twice as much on the horizon as in the for- 
mer case. 
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FIG. 6: Misner-Sharp, Brown- York, and Liu-Yau masses 
for a cluster of galaxies with NFW density profile in a flat 
background. The total density parameter Q = p/p cr as a 
function of the physical radius is also shown. Horizon radius 
is about O.Olpc. 
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FIG. 8: The Misner-Sharp and Brown- York matter flux com- 
puted along apparent horizon. Note the rapid increase of the 
matter flux rate up to the 'freezing point' point of the Horizon 
and the following decrease. 
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FIG. 9: The apparent horizon line in t-R plane. 



FIG. 7: Misner-Sharp, Brown- York and Liu-Yau masses for 
a galaxy with NFW density profile in a flat background. The 
density parameter f2 = p/p cr is also shown. Horizon radius is 
about I0~ 4 pc. 

Fig. (fT0|) compares the BY mass for two different NFW 
density profiles at a given time after the onset of the 
singularity. It shows the increase of the horizon and black 
hole or structure masses, and their shift to the right with 
the increasing density profile. 



V. DISCUSSION AND CONCLUSIONS 

Being faced with the challenge to define mass of struc- 
tures in a dynamical setting within general relativity, 
we have used models of mass condensation within a 
dynamical cosmological background to gain concrete 
insights of the similarities and differences between some 



of the more familiar mass definitions. The first two toy 
models representing collapsed overdense regions within 
asymptotically open-flat and closed-flat FRW models, 
show similar behavior at distances not very far from 
the collapsed region. The Brown- York mass, however, 
becomes negative for the open-flat model (FigllJ at large 
distances relative to the place of the void where the 
density reaches the background one. For both models 
the BY mass tends to zero at infinity, as expected. It 
seems that independent of any density profile as the 
initial condition, if we wait enough, there is always a 
void before the density reaches the background value. 
At about the same distance, there is always a maximum 
of the mass, the 'structure mass maximum', which in- 
creases and move to the larger distances from the center 
as the density profile increases with the time through 
more infall of matter to the singularity. At the central 
singularity for R — the Brown- York mass is zero 
in accordance with the BY mass of the Schwarzschild 



9 



10 100 1000 10000 100000 1000000 

' l.E+04 

Mbryoar ■ 'Omega omega 



^ l.E+14 - 
I 





0.1 1 10 100 1000 10000 100000 1000000 

R(Mpc) 

FIG. 10: Behavior of Brown- York mass for two different 
density profiles at the same time. 



metric but in contrast to the Misner-Sharp, Hawking, 
and Hay ward mass. At distances less than the void 
the BY mass is less than the Misner-Sharp mass. It 
is interesting to note that all mass definitions lead to 
almost the same value near the structure maximum. 
In the case of the more realistic NFW density profile 
these features are even more distinguished. As can be 
seen from the Fig. (J6j) , the structure maximum moves 
to the left at distances of the order of few Mpc for 
masses of the order of clusters. This is the same physical 
distance at which all three masses almost coincide and 
are equal to each other. At distances above lOMpc 
the LY and MS mass are almost equal except for a 
maximum LY mass at distances corresponding to the 
cosmological BY maximum mass and differing from the 
Misner-Sharp ones. The BY mass at distances larger 



than place of the void may be defined as a function of 
the comoving radius or physical radius, being almost 
the same. It differs, however, substantially from the LY 
and Misner-Sharp mass, and remains almost constant 
up to large cosmological distances before reaching the 
last cosmological maximum. 



The results obtained so far is indicative enough that 
the mass definitions may differ substantially. It is not 
said, however, that it may have any impact on our as- 
trophysical mass determinations. In fact the mass defi- 
nition in astrophysics is not as trivial as it may seems in 
the Newtonian dynamics. It is not even clear that there 
is any need at all for the concept of 'mass' in any astro- 
physical or cosmological setting. We are currently using 
our models to see if one can see any discrepancy between 
the general relativistic mass definitions and the one in 
Newtonian dynamics. Given the the Newtonian approx- 
imation for weak fields in general relativity, it is a legiti- 
mate question if this limit may also be used at cosmolog- 
ical distances (29[. We have already seen how different 
Misner-Sharp mass is relative to BY one! To tackle such 
questions we are currently applying our model to grav- 
itational lensing phenomena Q and the rotation curves 
of point masses within dynamical structures to see any 
deviations from the Newtonian approximation. Specifi- 
cally, we model a lens as a structure in the cosmological 
background and solve the geodesic equations numerically 
in a general relativistic framework using our model struc- 
ture. Note that the structure maximum mass occurs at 
points where the density is of the order of 10~ 29 g/cm 3 
and the gravity is weak enough to assume the Newto- 
nian approximation. It is, however, not clear that we 
can ignore the nonlinear effects of general relativity at 
such large distances [29j . 
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